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Fundamental aspects of ultracold collisions between identical bosonic or fermionic dipoles are 
studied under quasi-two-dimensional (Q2D) confinement. In the strongly dipolar regime, bosonic 
and fermion species are found to share important collisional properties as a result of the confining 
geometry, which suppresses the inelastic rates irrespective of the quantum statistics obeyed. A 
potential negative is that the confinement causes dipole-dipole resonances to be extremely narrow, 
which could make it difficult to explore Q2D dipolar gases with tunable interactions. Such properties 
are shown to be universal, and a simple WKB model reproduces most of our numerical results. 
In order to shed light on the many-body behavior of dipolar gases in Q2D we have analyzed the 
scattering amplitude and developed an energy-analytic form of the pseudopotentials for dipoles. For 
specific values of the dipolar interaction, the pseudopotential coefficient can be tuned to arbitrarily 
large values, indicating the possibility of realizing Q2D dipolar gases with tunable interactions. 
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Extensive experimental and theoretical efforts have re- 
cently been devoted to explore the production of dipo- 
lar gases and to uncover novel quantum phases resulting 
from the anisotropic nature and long-range character of 
the dipolar interaction [l[ . The recent experimental real- 
ization of a dense sample of ultracold ground-state KRb 
molecules Q has opened up ways to realize several new 
phenomena, ranging from ultracold chemistry, condensed 
matter physics and quantum information [2hJ] . A crucial 
step towards realizing such ideas is to understand the col- 
lisional properties of dipoles at ultralow temperatures. In 
fact, recent theoretical work has proposed a classification 
of ground-state molecules in terms of their chemical re- 
activity [f| . A number of such molecules are expected to 
be extremely reactive while others have energetically for- 
bidden decay channels. Although reactive molecules are 
great candidates for studying chemical dynamics 0, 0] , 
such processes limit the lifetime of the gas. The introduc- 
tion of a quasi-two-dimensional confinement (Q2D) adds 
a new handle for controlling the physics of the problem 
0] , which drastically changes this scenario and ul- 
timately leads to improved stability. This should in turn 
enable an exploration of various many-body phenomena. 

In this paper, we explore the effect of Q2D confinement 
on collisions of bosonic and fermionic dipolar species. We 
have found that in the strong dipolar regime the ultra- 
cold scattering properties of Q2D dipoles do not depend 
on the details of the short-range interactions. The ef- 
fective interactions are also characterized by calculating 
the Q2D pseudopotential coupling constants. Our anal- 
ysis shows that at any finite collision energy there exists 
a value of the dipolar interaction at which the coupling 
constant diverges, a phenomenon that can be traced back 
to the Ramsauer-Townsend effect. Although our present 
analysis does not address the issues of implementing our 
pseudopodential formalism to the many-body problem, 
we believe that our findings can shed light on the possi- 



ble ways of controlling the interactions in a dipolar gas. 

In order to study dipolar collisions in Q2D, an adia- 
batic separation of the radial and angular motion is in- 
troduced. This expresses the wavefunction as ty(r) = 
e i™$Y! ll ,r~ 1 F v (r)$v(r',0), where r, 0, and (f>, are the 
spherical radius and angles, respectively, and m is the 
angular momentum projection. Here, v is the channel 
index, F v is the z/-th radial wave function and $ is the 
channel function. In the adiabatic representation, the 
Schrodinger equation reduces to a simple system of ordi- 
nary differential equations given (in atomic units) by, 
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F„+Y, W m ,{r)F v , =0, (1) 



where fi is the two-body reduced mass, E is the total 
energy, W vv > = - (P vv >d/dr + Q vv >) /2fi, with P vv , = 
($„|3/0r|$„) and Q vu , = {§ u \d % /dr 2 \$ v ,), are the 
nonadiabatic couplings which drive inelastic transitions, 
and W v = U u — Q uv /2ii are the adiabatic potentials sup- 
porting bound and quasi-bound states. The adiabatic 
potentials and channel functions are obtained by solving 
eigenvalue equation for the angular motion: 
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solved for fixed values of r and in concert with the proper 
bosonic and fermionic boundary conditions at 9 = tt/2. 
In the equation above, v sr (r) = Dsech 2 (r/r ) is the 
short-range isotropic potential, where r is the charac- 
teristic range of the interactions (similar to the van der 
Waals length) , and the next two terms are the anisotropic 
dipole-dipole and Q2D confinement potentials, respec- 
tively. The dipolar interaction and confinement intro- 
duce two new important length scales into the system, 
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namely, the "dipolc length" di = [i.d 2 m , where d m is the 
dipolc moment in a.u., and the reduced- mass oscillator 
length a\ Q = fiWho, where Who is the harmonic trap fre- 
quency experienced by each molecule in the z direction. 

In Figure [T] we show a set of adiabatic potentials ob- 
tained by solving Eq. for the lowest bosonic (m = 0) 
and fermionic (m = 1) symmetries with even z-reflection 
parity. For distances r < a/j , the adiabatic potentials 
are mainly controlled by the short-range interactions (v sr 
and dipolar). As r increases the confinement becomes im- 
portant and the transition between 3D and 2D physics 
is characterized by the series of avoided crossings along 
the values of the averaged confinement potential r 2 / £ta^ 
(see the solid curve in Fig. [T]). For ultracold collisions, 
however, the scattering properties are mainly determined 
by the asymptotic forms of the adiabatic potentials which 
we found, for r 3> r; lo = (Aa^dt) 1 ' 5 , to be given by 

v 2 -1/4 di 
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Here E^J = (n Iy + l/2)//xo 2 lo is the oscillator energy level 
and n v is the principal quantum number. Notice that in 
contrast to the 3D case [9| the dipolar interaction at large 
r is repulsive 1/r 3 for both bosonic and fermionic species. 
This fact has important consequences for the collisional 
properties of dipoles in Q2D, as we discuss below. 
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FIG. 1: (color online). Adiabatic potentials for Q2D bosonic 
(solid) and fermionic (dashed) dipoles (see text). 

Having determined the adiabatic potentials and cou- 
plings, we calculate elastic and inelastic rates by solv- 
ing Eq. ([T]) and matching our numerical solutions to the 
proper cylindrical asymptotic solutions to obtain the cor- 
responding T-matrix and the collision rate coefficients Q , 

. , K in = ^£A m |T^| 2 ,(4) 

where A m = 2 — 8 m $ and i and / label the initial and 
final collision channels. Here, we study collisions from the 
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lowest transverse mode (ni = in Fig. [T]) and we have 
added an artificial deeply bound channel (with adjustable 
coupling) to incorporate the effects of inelastic collisions 
relevant for reactive ground-state polar molecules 

In general, a calculation of elastic and inelastic rates 
in the strong dipolar regime, kdi 3> 1, where k 2 = 
2(i{E — ErJ), requires the inclusion of a large number of 
partial waves in Eq. Q. Here, however, we have found 
that the rates in this regime assume a universal form, al- 
lowing for a simple WKB model for the phase-shifts that 
tremendously economizes the calculations. We compare 
our numerical results for the elastic rates for m = and 
1 with the WKB model in order to demonstrate this uni- 
versality and present the calculation for the total rates 
[Eq. 0] for values of m up to 200 using our WKB result. 
The topology of the adiabatic potentials (Fig. [l| allows 
us to directly apply the WKB phase-shift formula derived 
in Ref. [l(J. Here, however, we write the phase-shift as 
a sum of two terms, 
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where 5 l £ = [K m (r) — k]dr + — fcr^ is the long- 
range phase-shift, r c m being the classical turning point 
111, and 6% = tan _1 [ie _: 
phase-shift. Here K m (r) 



tan <F 

\k 2 - 



is the short-range 
(2di + m 2 r)/r 3 ]V2, 

exp[— 2 J r ™ \K m (r)\dr]. In our formulation, the 
only non-universal part of the phase-shift comes from 5^ 
through the phase accumulated for r < r^ D . We have 
determined, however, that for kdi 1, the phase 5^ is 
exponentially small [j m oc (df/a^o) 2 / 5 ], due to the 1/r 3 
barrier in the entrance channel [Eq. (|3J)], except in an 
extremely narrow region near a dipole-dipole resonance 
[l2| . As a result the scattering problem becomes univer- 
sal, depending only on the long-range physics encapsu- 
lated in 5^ = 5 l ^(k,di). The fact that, in Q2D, dipole- 
dipole resonances are extremely narrow, seems likely to 
make the exploration of Q2D dipolar gases with tunable 
interactions prohibitively difficult. However, as we will 
show later, tunability is still possible, although instead 
due to the physics of the Ramsauer-Townsend effect. 

Figure [2] demonstrates some of these points through 
our numerical calculations (for m = and 1) and WKB 
results for the partial rates KJ^ = (8ft//i)A m sin 2 S m (k) 
and /C™ for a/ lo /r = 5 and k = 0.071/r (vertical lines in 
Fig. [2]). Figs. [2] (a) and (d) compare our numerical results 
for m — and m = 1 (thick dashed lines) with the WKB 
results from Eq. ((5]) (thin dot-dashed lines). Extremely 
good agreement is found for kdi > 1, but for kde < 1, 
this agreement deteriorates and dipole-dipole resonances 
[l2T | becomes visible. In fact, throughout the entire range 
of values of di there exist dipole-dipole resonances, but as 
we mentioned above, they are extremely narrow and only 
become visible in the inelastic rates [Figs. [2] (b) and (e)]. 
Interestingly, all partial rates oscillate as di increases, as 
a result of the Ramsauer-Townsend effect. Evidently, 
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FIG. 2: (color online). Elastic and inelastic rates for bosonic, (a) and (b), and fermionic, (d) and (e), dipoles as a function of 
di for a,ho = 5ro and k = O.OTlr,^ 1 . In the strongly dipolar regime, kdt 3> 1, the elastic rates approaches to their semi-classical 
expectation 0] [solid straight line in (a) and (d)] and the inelastic rates are strongly suppressed [(b) and (e)]. (c) and (f) shows 
the scattering amplitude displaying anisotropy effects (see text). 



the total rates will not display actual zeros since the 
partial rates are out of phase. Nevertheless, the total 
rates, shown in Figs. [5] (a) and (d) as solid lines, oscillate 
around the semi-classical result of Ref. Q (solid straight 
line) . Still, such oscillations are universal and reflect the 
Ramsaucr-Townscnd effect. Figures [5] (b) and (e) also 
show another dramatic effect due to the Q2D confine- 
ment. The inelastic rates for both bosonic and fermionic 
dipoles are suppressed as exp[— 2(de/ai lo ) 2 ^ 5 ] for values 
of di beyond a^ , in agreement with Ref. 0. This result 
emphasizes the importance of the confinement in allow- 
ing for the realization of stable ultracold dipolar gases, 
irrespective of their quantum statistics. 

Figures [5] (c) and (f ) shows other properties of dipo- 
lar scattering through plots of the scattering amplitude 
versus the azimuthal angle, (f>, and di. As one can see, 
the dipolar scattering is by itself anisotropic. For both 
bosonic and fermionic dipoles [Figs. [2] (c) and (d)], the 
preferential outgoing scattering flux occurs around </> = 0, 
7r, and 2n. Therefore, two dipoles tend to scatter in such 
a way that they can either completely repel or else sim- 
ply "ignore" each other. Nevertheless, our results also 
show that there exist other specific directions in which 
dipoles can scatter. These are indicated by the less in- 
tense "fringes" between the maximum values of the scat- 
tering amplitude. From a more general perspective, the 
properties of the scattering amplitude could prove useful 
in many-body theories of dipolar gases fl3| . 

Another way to analyze the many-body behavior is 
through the pseudopotentials for Q2D dipoles 14, [ljj. 
However, pseudopotentials for dipoles (both in 3D and 
Q2D) suffer from a number of still-unresolved prob- 
lems in the limit k — > 0, and thus their applicabil- 
ity in many-body contexts has remained questionable. 
Here, we propose a modified version of Q2D potentials 
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FIG. 3: (color online), (a) Q2D coupling constant, g m , and 
(b) its equivalent energy-analytic form, g^- Figure shows 
interesting "resonant" behavior related to the Ramsauer- 
Townsend (see text). 



which eliminates the non-analyticity of the pseudopoten- 
tials in the limit k — 5> 0. Our approach consists in de- 
termining the pseudopotentials with respect to the en- 
ergy analytic radial 2D solutions in the spirit of gener- 
alized quantum defect theory In the present case 
this amounts to replacing the usual regular and irregu- 
lar Bessel solutions by a new pair that are both entire 
analytic functions of energy: = (kb)~ m J m (kr) and 
G a m = (kb) m [Y rn (kr) - l\n(kb)J m {kr)l where J m and 
Y m are the usual Bessel functions and b is an appropriate 
length scale which, based on our numerical calculations, 
is set to b = d(,. From this analysis, and following 14J, 
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we obtain the energy-analytic pseudopotentials as, 
«»(r) g a m (k) x S(r)R a m (r), 



(6) 



where 3^ is the energy analytic coupling constant and 
i?^ is the regularization operator given, respectively, by 
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In the equations above, /3 m = e~ 7+ ^p =1 p , where 7 ss 
0.577 is the Euler constant and <5£, is the energy analytic 
phase-shift, which relates to the physical phase-shift by 



tanC(fc) = 



tan S m (k)/(kdt) 



l-l{kd l ) 2m Hkd l )t & rx5 m {k) 



(10) 



Note that our form of the pseudopotcntial has an ex- 
plicit dependence on the dipole length dg while most of 
the k dependence comes through tan<5^. In fact, for 
m = our pseudopotential becomes energy independent 
in the kdg <C 1 regime. These properties appear to be ab- 
sent from previous analyses of Q2D pseudopotentials but 
they appear to be desirable for many-body calculations. 
In Figure[3]we show the contrast between using the usual 
form of the pseudopotentials and the energy-analytic 
form [Eqs. ©-©]• Figures [3] (a) and (b) shows our re- 
sults for the physical coupling constant g m [obtained by 
replacing dg by 1/fc and 8^ by S m in Eq. Q] for the 
same set of parameters chosen in Fig. [2] As one can see, 
g m increases for large m t&n5 m /k 2m oc l/fc 2 " 1-1 (m 7^ 0), 
while [Figs. [3] (c) and (d)] decreases for kdg ^> 1. This 
indicates that our form of the pseudopotcntial might be 
more suitable in the strongly dipolar regime. Apart from 
the above considerations, we note that for every m, g m 
diverges at particular values of dg [l5| . These divergences 
are related to the Ramsauer-Townsend effect and do not 
correspond to a zero-energy bound state. If these "res- 
onances" persist in a many-body context it can indicate 
that the Q2D dipolar gas could undergo a transition to 
a strongly correlated regime where the interactions can 
be tuned from repulsive (g m > 0) to attractive {g m < 0) 
almost at will by varying the external electric field. More- 
over, the positions of the poles in g m , according to our 
model, should be universal in the regime kdg 3> 1. How- 
ever, since the positions of resonant features depend on k 
they might be smeared out in a quantum gas. But they 
could be observable in a collision between two clouds at 
a definite relative momentum [T3|. Now, the fact that 
the poles in g^ are visible only for m = emphasizes 
the necessity of extending the many-body analysis in the 



energy-analytic framework, a task beyond the scope of 
our present investigation. 

In summary, we have studied scattering properties of 
bosonic and fermionic dipolcs under Q2D confinement. 
Both species display universal behavior, which signifi- 
cantly improves the prospects of creating a stable gas 
of Q2D dipoles. We have also computed two-body pa- 
rameters, such as the scattering amplitude and coupling 
constants, which can be important for many-body treat- 
ments of strongly interacting dipolar gases. 
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